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ABSTRACT 


Vertical  distributions  of  suspended  matter  which  has  been  observed 
by  Ewing  and  Thorndike  at  three  deep  stations  on  the  continental  sLpe  south¬ 
east  off  New  York  with  a  new  nephelometer  are  discussed  frcm  a  point  cf 
view  of  turbulent  diffusion  of  sediment  particles.  It  is  shown  that  theoreti¬ 
cal  curves  obtained  with  eddy  diffusivity  which  vanishes  at  the  bottom  and 
the  top  of  a  nephelcid  layer,  varying  as  a  quadratic  function  of  depth,  do 
not  agree  with  the  observed  distributions,  though  the  model  of  eddy  diffusi¬ 
vity  has  been  proved  to  be  suitable  for  suspended  particles  in  an  open 
channel.  It  is  found  that  eddy  diffusivity  proposed  by  Rossby  and  Montgomery 
for  the  lower  layer  of  atmosphere  yields  vertical  distributions  of  concentra¬ 
tion  similar  to  those  observed.  The  eddy  diffusivity  determined  from  curve 
fitting  techniques  consist  of  three  parts:  a  lowest  part  linearly  increasing 
with  distance  from  the  bottom,  a  middle  part  quadratically  decreasing  with 
distance  &nd  an  upper  part  of  a  constant  value.  When  settling  velocity  of 
suspended  matter  is  taken  as  3  x  10  (cm/  sec),  the  ranges  of  eddy  diffusi¬ 
vity  determined  at  three  stations  are  0. 1  to  1  (cm^/sec)  for  the  upper  part 
and  0,  6  to  36  (cm^/ sec)  for  the  maximum  value  at  the  top  of  the  lower  part. 
Speculation  on  causes  of  transparent  zones  is  presented.  Origins  of  turbu¬ 
lence  which  produces  boundary  layers  in  adjacent  to  the  bottom  and  residual 
eddy  diffusivity  in  the  upper  part  of  the  nepheloid  layer  are  discussed. 
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lw  INTRODUCTION 


Recent  efforts  to  measure  currents  in  the  great  depths  of  the  ocean  have 
shown  this  to  be  one  of  the  most  difficult  tasks  in  oceanography.  A  more 
formidable  task  still  is  the  measurement  of  turbulence  structure  near  the 
deep  sea  bottom.  In  fact,  it  is  more  difficult  to  set  up  current  measure¬ 
ment  devices  in  the  deep  water  than  in  the  shallow  water  in  order  to  record 
fluctuations  of  currents  for  calculating  the  energy  spectrum  of  turbulence. 
Also,  there  is  no  reliable  instrument  for  recording  fluctuations  of  tempera¬ 
ture  or  salinity  in  the  deep  sea  owing  both  to  engineering  difficulties  and  to 
a  condition  that  gradients  and  fluctuations  of  these  quantities  are  so  small 
in  the  great  depth.  Lack  of  gradients  of  temperature,  salinity  and  concentra¬ 
tion  of  other  chemical  components  in  the  deep  sea  reduces  the  validity  of 
indirect  methods  of  determination  of  eddy  diffusivity  from  distributions  of 
these  quantities,  because  in  such  methods  not  only  the  gradients,  but  also 
the  second  order  derivatives  of  the  concentration  with  space  coordinates 
are  necessary.  (Ichiye,  1962)  It  is  almost  out  of  the  question  in  the  present 
state  of  the  art  to  use  dye  diffusion  techniques  which  have  been  utilized 
successfully  in  the  shallow  water  and  in  the  upper  layer  of  the  deep  sea, 
considering  the  immense  difficulties  in  release  and  detection  of  small  dye 
patches  in  the  great  depth. 

On  the  other  hand,  importance  of  information  of  turbulence  in  the  deep 
ocean  has  been  realized  in  oceanography  in  recent  years.  Particularly, 
knowledge  of  eddy  diffusion  is  essential  to  disposal  of  waste  including 
radioactive  materials  in  the  deep  sea  and  dispersion  of  fallout  from  nuclear 
testing.  So  far,  only  three  authors  have  estimated  eddy  diffusivity  in  the 
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deep  part  of  the  Philippine  Trench  and  South  Atlantic  Ocean  by  use  of  indirect 
methods  whose  validity  is  doubtful  as  stated  before  (Defant,  I960). 

Ewing  and  Thorndike  (1964)  have  recently  determined  vertical  distribu¬ 
tions  of  suspended  matter  near  the  bottom  of  the  continental  slope  off  the 
Atlantic  Coast  of  the  United  States  by  use  of  a  new  nephelometer.  These 
distributions  together  with  other  data  yield  rather  accurate  estimations  of 
eddy  diffusivity,  Since  the  instrument  can  be  operated  to  any  depth  of  the 
ocean,  it  will  give  a  simple  method  to  determine  eddy  diffusivity  in  the 
great  depth. 

2.  CALCULATION  OF  EDDY  DIFFUSIVITY 

In  rivers  and  laboratory  flumes,  vertical  distributions  of  suspended 

sediment  have  been  utilized  to  determine  turbulence  structure  of  the  flow. 

O'Brien  (1935)  and  Hurst  (1929)  considered  the  equation  for  a  concentration 

C  for  a  constant  eddy  diffusivity  A. 

w  C  +  A  _dC_  =  0  (1) 

dz 

where  w  is  settling  velocity  of  the  sediment.  Equation  (1)  can  be  inte¬ 
grated  at  once  to  give  the  concentration 
C  =  exp 

in  which  Ca  is  the  concentration  at  some  arbitrary  level  z  =  a  , 

Rouse  (1937)  assumed  the  similarity  between  the  eddy  diffusivity  A  and 
eddy  viscosity  V  ,  the  latter  of  which  is  given  by 


1-* 


-1 


(  z  -  a) 


i 


(2) 


V  skcu^z(l-z/h) 


(3) 
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when  it  is  assumed  that  the  vertical  profile  of  the  current  velocity  obeys 

Prandtl's  logarithmic  law.  In  Equation  (3)  ko  is  Harman' s  constant,  h 

is  the  depth  of  the  flume  and  u*  is  the  shearing  velocity  defined  by 

u*  =  kQ  z  (du  /  dz )  (4) 

where  u  is  the  velocity.  With  the  assumption  that 

A  =  p  v  (5) 

integration  of  Equation  (1)  yields 

C/C a  =  /  h  -  z  a  YV  (6) 

V  h  -  a  z  ) 

where 

y  =  w  (  p  k0  u*  )  ~1  (7) 

In  order  to  determine  the  eddy  diffusivity  from  the  vertical  distribution 
C  by  use  of  (6),  the  shearing  velocity  u*  must  be  obtained  from  (4)  by 
measuring  the  shear  and  the  coefficient  Y  must  be  determined  by  curve  fit¬ 
ting  method.  However,  it  is  obvious  that  the  distribution  of  suspended 
matter  in  the  deep  sea  does  not  correspond  to  the  model  which  yields 
Equation  (6)  as  indicated  in  Fig.  1,  in  which  typical  readings  of  the  Thorndike 
nephelometer  as  the  concentrations  are  plotted  against  the  distance  froxn  the 
bottom  besides  the  theoretical  curves  of  (6)  for  different  values  of  Y  .  The 
data  plotted  in  Fig.  1  were  obtained  at  Station  RC  8-5  which  is  located  at 
about  37.5°  N  and  70°  W  and  4023  meters  deep,  as  shown  in  Fig.  3  (Ewing 
and  Thorndike,  1964)  The  theoretical  curves  are  calculated  by  taking 
h  =  431  m  and  a  =  1  m,  considering  the  fact  that  the  reading  of  the  nephelometer 
showed  the  minimum  value  of  -0.2  throughout  above  431  m  from  the  bottom. 
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Since  it  is  unbelievable  that  the  logarithmic  law  for  the  velocity  profiles  is 
valid  throughout  this  entire  layer,  the  conspicuous  disagreement  between  the 
observed  curves  is  natural. 

The  observed  distribution  of  Fig.  1  may  be  obtained  from  the  theoretical 
relationships  (Z)  by  taking  two  different  values  of  A  which  can  be  determined 
by  fitting  two  straight  lines  to  the  observed  curve  plotted  in  semi-logarithmic 
scale  near  the  bottom  and  near  the  upper  limit  of  the  nepheloid  zone.  However, 
this  process  yields  the  discontinuous  distribution  of  vertical  eddy  diffusivity. 

On  the  other  hand,  Rossby  and  Montgomery  (1935)  developed  a  theory  on 
turbulent  shear  flow  in  the  lower  atmospheric  layer  based  on  Prandtl' s 
boundary  layer  theory.  According  to  them,  the  eddy  viscosity  increases 
linearly  with  height  from  the  surface  to  the  upper  limit  of  the  frictional  boundary 
layer,  yielding  the  logarithmic  profile  of  wind  speed  and  then  it  decreases  to 
zero  at  the  upper  limit  of  the  planetary  boundary  layer,  in  which  the  wind 
vector  changes  with  height  similar  to  Ekman's  spiral.  It  turned  out  that 
this  modal  of  the  eddy  diffusivity  also  does  not  yield  the  vertical  distribution 
of  the  concentration  as  observed  in  the  nepheloid  zone.  However,  if  the 
existence  of  a  layer  of  constant  eddy  diffusivity  is  aiesumjsd  in  iie  upyer  part  of  the 
turbidity  layer  as  susggested  by  points  scattered  on  a  straight  line  in  Fig.  1, 
the  theoretical  curve  can  well  be  fitted  to  the  observed  one.  This  method 
will  be  shown  next. 

The  eddy  diffusivity  A  is  considered  to  consist  of  three  parts: 


A  =  Aq  =  const 

(  z  >  D  ) 

(8a) 

A  =  kj  (  h  -  z  )2 

(  D  >  z>  H  ) 

(8b) 

A  =  k2  z 

(  H  >  z  >  0  ) 

(8c) 

: 


. 


-5 


The  distributions  expressed  by  (8b)  and  (8c)  are  the  same  as  proposed 
by  Rossby  and  Montgomery  (1935).  The  relation  (8a)  is  introduced  by  con¬ 
sidering  the  observed  curve  as  stated  above.  The  coefficients  kj  and  k2  can 
be  determined  by  continuity  of  A  at  z  =  H  and  z  =  D.  The  distribution  of 
concentration  C  is  obtained  by  integrating  Equation  (1)  with  z  for  values  of 
A  given  by  (8a),  (8b)  and  (8c),  When  the  distance  from  the  bottom  is  expressed 
by  a  non-dimensional  quantity  s  =  z  /  D,  the  distribution  of  C  is  given  by 


C 

C 

C 

whe  re 


C_,  exp  )  -  m(s  -  1  )  l  (  for  s  >  1  ) 

i  > 

cd  exP  ~  m(sh  -  l)  “  (  sh  “  1)Z  (sh  “  j 

ca  (sa  /  6  )  a  for(sh>s  >  sa  ) 


(9a) 

for  1  >  s  >  s^  (9b) 

(9c) 


m 


=  w  D  /  A  ,  sh  =  h  /  D, 


H 


=  H/D,  sa  =  a/D, 


a  =  sH  (  sh  -  1)' 


x  (sh  -  sH  ) 

and  and  Ca  are  the  values  of  C  at  z  =  D  and  z  =  a,  respectively.  The 
eddy  diffusivity  of  Equation  (  8b)  and  (8c)  can  be  expressed  by 

A  =  Aq  (  sh  -  s  )2  (sh  -  1  )-2  (for  1  >  s  >  sH  )  (8b1) 

A  =  A0  sH  s  (  sh  -  sH  )  (8c1) 

The  numerical  constant  w  /  AQ  can  be  determined  from  the  slope  C 
represented  in  semi -logarithmic  scale  as  in  Fig,  1  against  the  depth  for  the 
upper  linear  part  of  the  curve.  The  value  of  D  is  determined  as  the  lower 
depth  of  the  linear  part  of  the  curve  which  is  represented  by  Equation  (9a), 
The  reference  distance  a  is  arbitrary,  and  is  taken  as  1  meter,  which  is 
not  very  far  from  a  truth,  considering  some  observations  indicating  existence 
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of  slime  with  a  few  feet  depth.  Two  unknown  depths  h  and  H  can  be  deter¬ 
mined  by  fitting  the  curves  computed  from  (9b)  and  (9c)  with  the  observed  one. 
The  numerical  values  of  these  constants  thus  determined  are  as  follows: 

D  =  350  (m),  w /  Aq  =  2,  54  x  10“^  (cm” ^) 

Cd  =  0.36  (NU  or  Nephelometer  Unit),  Ca  =  1.17  (NU) 

sh  =  h  /  D  =  1,34,  sH  =  H  /  D  =  0.57 

The  value  of  corresponds  to  200  m  of  the  thickness  of  the  frictional  boundary 

layer.  The  theoretical  curve  of  C  computed  from  (9a)  (9b)  (9c)  for  these 
numerical  constants  is  plotted  in  Fig.  1.  It  is  seen  that  the  computed  curve  is 
in  good  agreement  with  the  observed  data.  It  may  be  argued  that  the  agreement 
is  obvious  because  several  numerical  constants  are  determined  by  curve  fitting. 
However,  it  should  be  remembered  that  there  are  a  number  of  degrees  of  free¬ 
dom  in  the  curve  fitting  and  particularly  the  agreement  in  both  curves  in  the 
regions  expressed  by  (9b)  and  (9c)  suggests  that  the  mechanism  of  suspension 
might  actually  correspond  to  the  proposed  theory.  The  vertical  distribution  of 
eddy  diffusivity  A  corresponding  to  the  theoretical  curve  shown  in  Fig.  1  is 
plotted  as  a  ratios  to  AQ  in  Fig.  2. 

In  order  to  determine  eddy  diffusivity  from  the  estimated  value  of  w  /  AQ  , 
we  have  to  know  the  settling  velocity  w  ,  which  can  be  computed  by  Stokes'  law, 
if  the  mean  diameter  of  the  suspended  matter  is  known.  On  the  other  hand, 
according  to  Arrhenius  (1963),  the  settling  velocity  of  the  suspended  matter  in 
the  ocean  is  several  orders  of  magnitude  larger  than  the  value  obtained  from 
Stokes'  law  for  particles  whose  diameters  are  of  the  order  of  l ^  or  less.  The 
settling  velocity  of  the  finer  fraction  of  the  suspenoid  (0.01  to  0.  5jJ  )  whose 
average  diameter  is  assumed  to  be  0.05  is  estimated  by  him  as  3  x  10“^ 
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(cm/sec).  Using  this  value  for  w  ,  we  obtain  A0  =  0.12  (cm^/sec)  with 
the  value  of  w  /  AQ  determined  from  the  observed  data  in  Fig.  1.  This  value 
seems  to  be  small  compared  with  values  of  eddy  diffusivity  estimated  by  others 
in  the  deep  sea  using  temperature  or  salinity  data  (Defant,  1961).  However, 
it  is  still  several  orders  of  magnitude  larger  than  the  molecular  diffusivity  of 
2  x  10**^  (cm^/sec).  Since  in  this  calculation  settling  velocity  of  the  finer 
fraction  of  particles  is  assumed,  the  actual  settling  velocity  may  be  larger  than 
the  assumed  one,  if  the  suspendoid  consists  of  mineral  particles.  All  this 
estimation  of  the  eddy  diffusivity  depends  on  accurate  value  of  w  .  Also,  it 
is  more  important  to  note  that  the  eddy  diffusivity  determined  from  the 
nepheloid  distribution  shows  maximum  at  about  200  m  from  the  bottom  than  to 
determine  numerical  values  of  the  eddy  diffusivity. 

The  vertical  distributions  of  A/Aq  for  Stations  RC  8-4  and  RC  8-6 
(Fig.  3)  are  also  determined  from  the  nephelometer  data  in  the  same  way  as 
for  Station  RC  8-5.  The  vertical  distributions  of  the  nephelometer  readings 
are  plotted  in  Fig.  4  with  the  theoretical  distributions  obtained  through  curve¬ 
fitting  processes.  For  Station  RC  8-6,  two  readings  taken  while  lowering  and 
raising  the  instrument  are  shown.  At  these  stations,  existence  of  transparent 
zones  makes  the  data  points  scattered  more  wildly  than  at  RC  8-6.  However, 
if  we  discard  the  transparent  zones,  agreement  between  the  theoretical  curves 
and  observed  values  is  good.  In  Fig.  5,  the  vertical  profiles  of  A/Aq  is 
plotted.  The  characteristic  quantities  for  theoretical  distributions  of  concentra¬ 
tion  and  eddy  diffusivity  are  as  follows: 

D  (m)  H  (m)  w  /  AQ  (cm-1)  AH  /  AQ 

RC  8-4  350  125  4.4  xl0“5  53  -  5 

RC  8-6  100  50  4.8  x  10**5 
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Apj  is  the  maximum  of  A  which  is  reached  at  z  =  H,  Comparison  of  the 
numbers  of  the  last  two  columns  in  the  above  table  with  those  of  RC  8-5  indi¬ 
cates  that  the  eddy  diffusivity  AQ  is  one  order  of  magnitude  larger  and  Apj 
is  one  to  two  orders  of  magnitude  larger  at  these  two  stations  than  at  RC  8-5. 

This  suggests  that  the  transparent  zones  which  are  abundant  at  these  two 
stations  are  related  to  intense  turbulence  in  the  nepheloid  layer, 

3.  Various  Factors  Pertinent  to  Estimation  of  Eddy  Diffusivity 

There  are  many  arguments  against  the  validity  of  the  Fickian  equation 
in  describing  the  diffusion  phenomenon  in  the  ocean  (Ichiye,  1962).  However, 
the  equation  may  be  valid  when  it  is  applied  to  the  process  which  occurs  in  a 
limited  area  for  a  steady  state  or  in  a  limited  period  of  time.  Also,  the  idea 
of  eddy  diffusivity  on  which  the  equation  is  based  is  not  inadequate  to  indicate 
the  state  of  turbulence  if  we  interpret  the  estimated  value  with  care.  Even  so, 
there  remain  many  problems  to  be  solved  in  order  to  determine  the  eddy 
diffusivity  by  use  of  distributions  of  suspended  matter. 

Firstly,  the  relationships  between  readings  of  the  nephelometer  and 
concentration  of  suspended  material  must  be  determined.  There  might  be  a 
saturation  concentration  beyond  which  the  nephelometer  gives  the  constant 
maximum  reading.  Secondly,  the  particle  sizes  of  suspended  matter  must  be 
measured  at  several  heights  from  the  bottom.  As  discussed,  in  the  end  of  the 
previous  section,  the  absolute  values  of  eddy  viscosity  cannot  be  determined  from 
the  distribution  of  concentration  without  knowledge  of  the  settling  velocity  of 
particles.  Also,  there  is  a  possibility  that  the  distribution  of  particle  sizes 
may  change  with  the  depth  and  thus,  the  concentration  merely  reflects  the 
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vertical  distribution  of  the  settling  velocity  with  a  constant  diffusivity. 

Even  if  the  validity  of  the  Fickian  equation  for  the  suspended  matter 
problem  is  established,  there  remains  the  doubt  about  applicability  of 
Equation  (1)  which  is  so  much  simplified  from  the  complete  Fickian  equation. 
Firstly,  the  assumption  of  steady  state  must  be  questioned.  Although  the 
nepheloroeter  readings  in  most  stations  show  the  same  values  within  an  accur¬ 
acy  of  the  instrument  while  lowering  and  hauling  up,  there  is  no  information 
on  the  change  of  concentration  with  time  longer  than  a  few  hours. 

Secondly,  neglect  of  the  convection  terms  due  to  water  movements  must 
be  checked.  It  can  be  shown  that  particularly  the  effect  of  vertical  velocity 
of  water  is  important.  The  Fickian  equation  in  a  steady  state  becomes 

u  8G/8x  +  v8C/9y  +  (wg  +  w)  8C/8z  =  Diffusion  terms  (10) 

where  u  and  v  are  components  of  horizontal  velocity  and  wg  and  w  are 
settling  velocity  of  particles  and  vertical  velocity  of  the  water,  respectively. 

So  far,  there  is  no  data  which  gives  exact  values  of  currents  and  horizontal 
gradient  of  the  concentration  and  thus,  we  have  to  be  satisfied  with  estimation 
of  orders  of  magnitude  from  the  existing  data.  The  order  of  magnitude  of 
horizontal  gradients  of  C  can  be  estimated  from  the  difference  of  vertically 
averaged  concentration  at  three  stations,  RC  8-6,  RC  8-5  and  RC  8-4  (Fig.  3) 
and  is  given  by 

Ord  (Hor.  Grad  C  )  =  0.1  NU/400  (km)  =  2.  5  x  10-9  (NU/cm)  <• 
where  NU  means  a  unit  of  the  nephelometer ,  The  averaged  vertical  gradient 
9  C  /  9z  is  estimated  from  the  vertical  distribution  of  C  at  one  station,  for 
instance,  at  RC  8-5  and  is  equal  to  1  (NU)  /  500  (m)  =  2  x  10**^  (NU  /  cm). 
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Therefore,  in  order  to  yield  the  same  effect  on  the  distribution  of  the  concentra¬ 
tion,  the  horizontal  component  of  the  velocity  must  be  10^  times  the  vertical 
component.  It  is  instructive  to  compare  this  figure  with  the  average  bottom 
slope  of  the  strip  of  400  km  width,  at  the  center  of  which  the  station  RC  8-5 
is  located.  The  slope  is  6  x  10**^  and  thus  for  the  current  flowing  parallel  to 
the  bottom  slope,  the  term  w8G  /  3z  is  one  to  two  orders  of  magnitude  larger 
than  the  first  two  terms  on  the  r.  h,  s,  of  Equation  (10), 

4,  Discussions 

Although  it  is  highly  speculative,  it  is  worthwhile  to  consider  causes  of 
transparent  zones  in  the  nepheloid  layer  at  some  stations  (RC  8-6,  for  example) 
reported  by  Ewing  and  Thorndike  (1964),  Since  the  effect  of  vertical  component 
of  the  water  movement  seems  to  be  predominant  on  the  distribution  of  suspended 
matter,  there  must  be  some  kind  of  vertical  motion  in  order  to  generate  the 
transparent  zones.  One  speculation  for  such  causes  is  internal  waves  of  the 
higher  modes  which  are  generated  due  to  the  stratification  caused  by  the  sus¬ 
pended  material.  Another  speculation  is  stationary  cellular  motion  of  convective 
origin.  The  vertical  streakiness  due  to  local  bottom  topography  of  the  current 
flowing  parallel  to  the  bottom  cannot  be  eliminated  as  a  cause,  since  such 
current  has  effective  vertical  components  as  explained  before.  If  data  of  the 
nepheloid  concentration  are  obtained  repeatedly  at  the  same  station  or  at  several 
stations  located  closely  to  each  other,  we  shall  be  able  to  determine  causes  of 
the  transparent  zones  as  well  as  the  details  of  bottom  current  from  the  change  of 
these  zones  with  time  or  space. 

The  vertical  distributions  of  eddy  diffusivity  determined  at  three  stations 
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indicate  the  presence  of  two  different  layers  (  z  >D  and  D>s>0  )  in  the 
deep  water  with  regard  to  turbulence  characteristics.  In  the  layer  between 
the  bottom  and  z  =  D,  eddy  diffusivity  varies  with  height  in  a  similar  way  to 
the  atmospheric  boundary  layer  near  the  ground.  This  similarity  seems  to 
be  superficial,  because  in  the  atmosphere  there  are  two  dominant  effects 
which  generate  turbulence:  shear  in  the  lower  layer  of  geostrophic  winds  and 
free  convection  due  to  differential  heating  on  the  ground  (Priestly,  1959)  while 
in  the  ocean  both  effects  seem  to  be  negligibly  small.  However,  a  quasi- 
geostrophic  nepheloid  current  which  will  be  discussed  elsewhere  (Ichiye,  1964) 
may  generate  turbulence  in  the  boundary  layer. 

Finite  eddy  viscosity  in  the  layer  above  D  was,  in  case  of  atmosphere, 
attributed  to  'residual  turbulence'  whose  energy  source  is  found  partly  in 
convective  motion  and  partly  in  shearing  forces  due  to  vertical  variations  in 
geostrophic  winds  (Rossby  and  Montgomery,  1935).  However,  the  convective 
motion  seems  to  be  very  small  in  the  great  depth  of  the  ocean.  Some  prelim¬ 
inary  experiments  on  effects  of  winds  and  gravity  waves  on  diffusion  of  dye  in 
water  were  made  in  a  racetrack-type  flume  of  23  cm  wide  and  50  cm  deep  of 
an  outward  periphery  of  a  major  axis  of  3.0  m  and  a  minor  axis  of  1, 5  m. 
(Ichiye,  1963)  .  The  results  indicated  that  the  dye  released  near  the  bottom 
was  diffused  much  more  rapidly  when  winds  were  blown  on  the  water  surface* 
than  when  gravity  waves  were  generated,  although  there  was  no  noticeable 
motion  near  the  bottom  in  the  former  case  compared  with  occurrence  of 
oscillatory  motions  in  the  latter  case.  This  suggests  that  wind  stresses  applied 
at  the  surface  may  become  an  effective  energy  source  of  turbulence  near  the 
bottom.  In  an  oceanic  scale,  energy  input  at  the  sea  surface  due  to 
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mete  or  o  logical  disturbances  is  mostly  consumed  in  the  upper  layer  above  a 
thermocline  and  only  a  small  portion  of  it  reaches  the  bottom.  However,  simpli¬ 
fied  analyses  for  a  two-layered  ocean  model  indicate  that  meteorological  dis¬ 
turbances  with  scales  and  duration  of  ordinary  fronts  and  cyclones  can  effectively 
generate  barotropic  currents  which  may  reach  the  bottom  (Veronis  and  Stommel, 
1956)  and  that  such  perturbations  in  the  deep  ocean  are  amplified  near  the 
western  coast  of  the  ocean  due  to  reflection  of  the  Rossby  waves  there  (Ichiye, 
1958).  It  is  necessary  to  measure  turbulence  in  the  deep  water  on  a  worldwide 
scale  in  order  to  determine  transport  of  turbulence  energy  from  meteorological 
disturbances  to  the  deep  water. 
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Distribution  of  nepheloid  concentration  (Nephelometer  Unit) 
versus  distance  from  the  bottom  (meter) 

Closed  circles:  Data  at  RC  8-5 

Broken  curves:  Theoretical  curves  computed  from  Equation  (6) 

Solid  curve:  Theoretical  curve  computed  from  Equations  (9a) 

(9b)  and  (9c).  A,  B  and  C  represent  portions 
of  the  curve  corresponding  to  (9c),  (9b)  and 
(9a)#  respectively. 


Vertical  profile  of  eddy  diffusivity  corresponding  to  the  theoretical 
distribution  of  nepheloid  expressed  with  the  solid  curve  in  Fig.  2. 
A,  B  and  C  correspond  to  each  portion  of  the  solid  curve  in 
Fig.  1. 


Locations  of  three  nephelometer  stations  (A)  and  the  bottom 
profile  of  the  transection  (B)  ,  (A  broken  line  in  (B)  indicates 
the  upper  limit  of  the  nepheloid  zone). 


Vertical  distributions  of  nepheloid  concentration.  (Closed  circles 
and  open  circles  indicate  the  readings  of  the  nephelometer  while 
lowering  and  raising  it,  respectively,  at  Station  RC  8-6.  Triangles 
indicate  the  readings  at  Station  RC  8-4  while  lowering  the  meter. 
Solid  and  broken  curves  indicate  theoretical  curves  obtained  from 
the  distribution  of  eddy  diffusivity  plotted  in  Fig.  5), 


Vertical  profiles  of  eddy  diffusivity  (A/Aq)  at  Station  RC  8-6 
(solid  curve)  and  at  Station  RC  8-4  (broken  curve). 
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